Searching for non-Fermi liquids under the holographic light 
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By expanding the set of background geometries beyond the commonly studied ones we identify 
those dual gravity models that may provide holographic descriptions for some prototypical non-Fermi 
liquid states of strongly correlated condensed matter systems. Specifically, we discuss prospective 
gravity duals of such iconic examples as the non-relativistic fermions coupled to gauge fields and 
Dirac fermions with the Coulomb interactions. 
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Quantum theory of strongly correlated fermions has 
long been in a rather desperate need for non-perturbative 
techniques, the use of which could allow one to pro- 
ceed beyond the customary (and often uncontrollable in 
the regime of interest) approximations when analyzing 
generic (non-integrable) systems. 

Despite all the effort, though, the overall progress has 
been quite limited. However, it's been argued that a re- 
cent proliferation of the ideas based on the conjecture 
of holographic duality^ may offer a possible way out of 
the stalemate, thus allowing the field to move towards 
a systematic classification of various 'strange' metallic 
(compressible) states that are commonly (and often in- 
discriminately) referred to as non- Fermi-liquids (NFLs). 

The holographic correspondence postulates a connec- 
tion between certain (non-Abelian, multi(iV)-component, 
and supersymmetric) field theory models and their grav- 
ity duals living in one extra dimension, so that for N 3> 1 
the strong-coupling regime of the former can be mapped 
on the weak-coupling one of the latter (and vice versa). 
The support for this general idea is provided by a host 
of circumstantial evidence gathered from the 'bona fide' 
theories of strings and hot QCD quark-gluon plasmas. 

Nevertheless, despite a gradually building confidence 
in the validity of the original holographic conjecture, the 
status of its recently proposed condensed matter appli- 
cations remains, by and large, unknown. 

Pursuing the phenomenological 'bottom up' approach, 
the initial studies produced a number of rather baffling 
results, which include multiple Fermi surfaces (merging 
into one critical 'Fermi ball' - or, rather, 'flat band'- in 
the extreme N — > oo limit), oscillatory frequency de- 
pendence and dispersionless poles of the fermion propa- 
gator, etc£. Some of those early findings have already 
been interpreted as spurious artefacts of taking the limit 
N — > oo, though. For one, no oscillatory dependence 
would arise in a more systematic 'top down' approach^ 
(besides, the same effect can also achieved by including 
a non- minimal (Pauli) fermion-gauge field coupling^). 

Putting aside the central question about a general ap- 
plicability of the original holographic conjecture to (typ- 
ically, neither non-Abelian/multicomponent, nor super- 
symmetric) condensed matter systems, one common lim- 
itation of the early studies was that their background 
metrics would typically be chosen from a handful of the 
well-known solutions to the classical Einstein-Maxwell 
equations. 



Obviously, such an 'under the light' search for holo- 
graphic NFLs lacks any physical input specific to a given 
strongly correlated system and, therefore, its chances of 
finding a gravity dual for that particular system appear 
to be rather hard to assess. 

In essence, the early studies have so far found just one 
type of the NFL behavior, dubbed 'semi-local criticality', 
whereby the fermion propagator features a non-analytical 
(and, in general, oscillatory) frequency dependence, but 
only a non-singular momentum one 

Here the function A(q) has simple zeros at each of 
the (potentially, multiple) 'Fermi momenta' q l F , whereas 
B(q) takes finite values at such points (see Refs.— for de- 
tails). 

This behavior suggests that at long distances the sys- 
tem effectively splits onto spatially uncorrelated 'quan- 
tum impurities', each of which exhibits a characteristic 
d = quantum-critical scaling. 

Although it was pointed out that the propagator (1) 
bears a certain resemblance to that expected in the con- 
text of some heavy fermion materials (see, e.g.,— and ref- 
erences therein), in the absence of any solid agreement 
with experiment that would suggest otherwise, this rem- 
iniscence might well turn out to be merely superficial. 

In view of such uncertainty, the task of putting the 
holographic correspondence on a firm ground and as- 
certaining the status of its predictions more definitively 
would be best achieved if this technique were applied to 
those situations where a prior insight has already been 
gained by some other means. To that end, it would be 
very helpful first to find the gravity duals of the already 
established (or, at least, suspected) NFL states. 

In what follows, we demonstrate that accomplishing 
this task requires one to extend the class of metrics well 
beyond the commonly studied classical examples. 

Prospective gravity duals of condensed matter systems 

The early applications of the holographic conjec- 
ture yielding Eq.(l) utilized the standard Reissner- 
Nordstrom (RN) 'black brane' solution which minimizes 
the Einstein- Maxwell action 
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Hereafter J = J dtdzd d xy / —detg fiL , stands for a covari- 
ant d + 2-dimensional volume integral, 1Z is the scalar 
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curvature, the second term represents a (negative) cos- 
mological constant of the asymptotically anti-dc-Sittcr 
space AdSd+2 of curvature radius L, k 2 is the Newtonian 
coupling, and e is the charge of the U(l) gauge field. 

The corresponding metric (throughout this paper, the 
speed of light c = 1) 



ds 2 = -f(z)dt 2 + g(z)dz 2 + h(z)dx 2 



(3) 



has non-zero components f(z)(z/L) 2 — (L/z) 2 /g(z) — 
1 - (1 + i M 2 )(z/z h ) d + 1 + fi 2 (z/z h ) 2d , h{z) = (L/z) 2 . Here 
fi is the (dimensionless) chemical potential of the fermion 
matter which is assumed not to react back on the metric, 
and Zh is the inverse radius of the horizon determined 
from f(z h ) = 0. 

The latter can remain finite even when the Hawking 
temperature T = (d + 1 — (d — l)n 2 )/4:TTZh vanishes, 
thereby giving rise to yet another artefact of the N — > oo 
limit, a seemingly non-vanishing entropy S(T — > 0) =/= 0. 
In the near-boundary z — > (ultraviolet or UV) regime, 
Eq.(3) recovers the standard AdSd+2 form, f(z) — g(z) — 
h(z) = (L/z) 2 . 

In fact, the 'locally critical' behavior (1) sets in at the 
extremal T — > limit 2 , where the near- horizon geometry 
approaches AdS2 x R d . Moreover, a similar (albeit more 
physically sound, entropy-wise: S(T) ~ T d ' v ) behavior 
was found for a variety of geometries which reduce to the 
one-parameter 'Lifshitz' metric f(z) = (L/z) 2 ' n , g(z) = 
h(z) — (L/z) 2 in the z 3> 1 (infrared or IR) regime^. It 
was also shown to result from an approximate (Thomas- 
Fermi) account of the fermions' back-reaction in the 
framework of the standard gravity (2), thus leading to 
the 'electron star' scenario— 

In fact, the latter metric naturally emerges alongside 
a whole class of more general solutions in the so-called 
dilaton gravity whose Lagrangian includes an additional 
bulk scalar field^ 
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In the minimal version, both, the dilaton potential 
U((f>) = d(d+ \)e 5<t '/L 2 and the effective gauge coupling 



Z(cf) 



— P i4> 



, are given by simple exponential functions 



with the coefficients S and 7. 

In what follows, we focus on the T = case and con- 
sider a still broader class of static and spherically sym- 
metric metrics 

f(z) = (L/zf\ g(z) = (L/z) 2a , h(z) = (L/zfP, 

(5) 

whereas at finite T one also has the freedom of altering 
the additional polynomial 'emblackening factor', similar 
to that in the RN solution (3). 

For any |3 / Eq.(5) can be reduced to a two- 
parameter family of metrics known as the 'hyperscaling 
violating' backgrounds 9 . The latter are characterized by 
the dynamical exponent £ and hyperscaling violation pa- 
rameter 9 

c = l±i^ 5 e = d^^- (6) 



which manifest themselves through the scaling properties 
of the excitation spectrum in the boundary theory: uj — > 
X^u for q — >• \q and that of the interval (3), ds —> \ e / D ds. 

In Refs. 9 , a number of 'top down' string scenaria was 
presented, by which the 'hyperscaling violating' geome- 
tries may arise. In that regard, the null energy criteria 
that must be obeyed by any physically sensible gravity 
dual with the metric (5) impose the inequalities 

/3(rj - 13 + a - 1) > 0, (rj - j9)(l - a + rj + dp) > (7) 

As in Refs.—, the generalized metrics (5) do not cover 
the z — > region and, therefore, require a proper UV 
completion. Thus, they should be viewed as gravity du- 
als of some effective IR field theories residing at a finite 
zq. Correspondingly, all the holographic propagators dis- 
cussed in the rest of this paper pertain to the renormal- 
ized operators from such effective theories, rather than 
those of the 'microscopic' boundary ones. The latter can 
be obtained from the former by virtue of the matching 
procedure, akin to that of Refs^. 
Semiclassical propagators 

In the holographic analyses the bulk fermions of mass 
m couple to the metric and gauge fields in the minimal 

way 

Sf = J $lti(idn + |[7A, lA^xv + eA v - m )V> (8) 

where 7 M are the 7-matrices and u^ v is the spin 
connection 1 . 

In the absence of explicit analytical solutions for the 
bulk fermion wavcfunctions in generic gravitational back- 
grounds, one can still resort to the semiclassical ap- 
proach. The equation for the Fourier-transformed wave- 
function ijj(z.uj,q) features an effective single-particle 
potential^ 



V(z) = m 2 + -j— + 



h(z) f(z) 



(9) 



which allows for two zero-energy solutions in the tunnel- 
ing region z < z < z t : 



l ( \ 1 ± J Zt dz^/g(z)V(z) 

i>±( z i UJ , ( i) ~ \n/^z) e (10) 

where the turning point z t is defined as V(z t ) = 0. 

Using Eq.(10) one then finds the effective IR theory's 
Green function as the reflection coefficient for the wave 
incident at z = z^ 



Gi R (u),q) = 



-4>+(z,u,q) 



where 



1 -a + P 



l-a + P 



S(LJ,q) = 2 dzy/g(z)V(z) 

J z 



(11) 
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Considering the metric (5) and focusing on the limit of a 
small fermion mass, one obtains the scaling behavior 



1 — a+7/ 



S(u,q)r.( — — ) 



(13) 



indicative of the underlying quasiparticle dispersion uj ~ 
governed by the dynamical exponent (6). 
In the complementary limit of a large mass, the semi- 
classical analysis can be more conveniently employed 
directly in the real spaced. In this regime, various 
quantum-mechanical amplitudes are dominated by the 
fermion paths that closely follow the classical trajectories 
(geodesies) obtained from the (imaginary-time) action 



S(t,x) = m / dzy/g(z) + f \z)(dr / 'dzf + h(z){dx/dz) 2 

(14) 

When evaluated upon such a trajectory, Eq.(14) reads 



S(t, x) 



dz\ 



/5(f) 
r(z) 



(15) 



where r(z) = 1 — U^/h(z t ) — V? T j f(z t ) is a function of 
the conjugate momenta H x = 5S/S(dx/dz) and IT r = 
5S/S(dr/dz) given by the integral equations 




dz g(z) 



f(z) y r{z) 



(16) 

and the turning point is obtained by solving the equation 
r(z t ) — 0. The minimal action (15) then controls the 
fermion propagator, G{t,x) ~ e~ s ( T ' x \ 

While an explicit computation of Eq.(15) can only be 
performed in some special cases, determining simpler, 
one-parameter, dependences S(t) and S(x) is possible 
for a broad variety of metrics. Specifically, for the metric 
(5) one obtains 



S(x)~ x *=^+e, S{t) 



(17) 



Notably, both Eqs.(12) and (15) elucidate the role of the 
radial variable z as an energy-like renormalization scale 
parameter. However, a direct correspondence between 
the two can not be readily established, since the Fourier 
transformation relating G(u>, q) and G(r, x) requires both 
functions, including their non-exponential prefactors, to 
be known across the entire ranges of their arguments. It 
might be possible, though, to relate their asymptotics by 
virtue of the saddle-point method, wherever applicable. 
Finite density fermions with singular interactions 
One important testing ground for the holographic con- 
jecture is provided by the theory of finite density fermions 
coupled to an Abelian gauge field. This problem has long 
been at the forefront of theoretical research where it was 
studied with a whole variety of techniques, although the 
case still remains unclosed. For instance, the recent re- 
sults of Refsii2 which revisited the attempts to obtain 



a self-consistent re-summation to all orders in the spirit 
of the Eliashberg theory^ indicate that a naive 1 /TV- 
expansion may not be as reliable as previously thought. 

Furthermore, long-ranged and retarded ('singular') in- 
teractions that allow for a similar description are often 
associated with the onset of ground state instabilities, 
and the concomitant NFL behaviors might occur even in 
those systems whose microscopic Hamiltonians involve 
only short-ranged couplings. 

Such interactions are mediated by gapless bosonic exci- 
tations of an emergent order parameter, and in all the di- 
verse reincarnations of the problem their gauge-like prop- 
agator conforms to the general expression 



D(u,q) 



1 



(18) 



Important pertinent examples include anomalous elec- 
tromagnetic skin effect in metals^, compressible Quan- 
tum Hall states with screened repulsive interactions^, 
critical spin fluctuations in itinerant ferromagnets^ and 
density fluctuations in 'quantum nematics' 17 , for all of 
which £ = l,p = 2. In contrast, normal skin effect 
and antiferromagnetic fluctuations would be described by 
£ = 0, p = 2, whereas compressible Quantum Hall states 
with the unscreened Coulomb interactions correspond to 
£=l,p=l. 

The asymptotic IR behavior of the propagator of 
fermions coupled to a gauge-like bosonic mode can be 
evaluated by means of the eikonal-type procedure^ 
which reduces the former to the phase factor taken along 
the classical trajectory 



G(t,x) ~< exp(i J Ap(z = z )dx t _ t ) > A -- 



(19) 



Here the averaging is performed over a (physical or ef- 
fective) gauge field governed by the propagator < 
A„Av >=D(w,q){6 tiv ~ q^qvj q 2 ), thereby resulting in 

. If diod d q . 1 — cos(wt — qx) . „. 

S ( T >*) = o / 7fW fl u ^ — r — W- 1 20 

2 J (2ir) a+L [iu) — vq)^ 

where v ~ qpii is the Fermi velocity in the direction of 
the vector x. 

To estimate the eikonal action (20) for a time-like in- 
terval 

, . If du;d d q . 1-coscjt P +i-d 
5(T) = 2 J (2^ D ^ q) (uJ^ - T ?+P (21) 

we first perform the momentum and then the frequency 
integrations, thereby discovering that at r —> oo the in- 
tegral is dominated by the frequencies ui ~ r _1 and mo- 
menta q ~ t _1 /'" +p . 

Moreover, the kinematics of fermion scattering is such 
that at small scattering angles (which is the regime 
amenable to the eikonal approximation) one finds u <C 
|vq| < |v x q| for all 2 - £ < d < 1 + p. Thus, the scat- 
tering momentum appears to be primarily directed along 
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the Fermi surface and perpendicular to the local Fermi 
velocity. 

With that observation in mind one can also compute 
the integral (20) for a space-like interval, thus obtaining 
the asymptotic behavior 



S(x) 



(22) 



It is then easy to see that the Eqs.(17) and (21,22) match, 
provided that the following relations hold 



1 



P- 



1-d 



1 — a 
1 — a + r\ 



P- 



d-1 
l-d 



P + Z 



(23) 



The above results pertain to the propagators of the ef- 
fective field theories which, unlike their underlying mi- 
croscopic counterparts, become essentially universal after 
having undergone renormalization down to the IR scale 
z . As such, they need to be contrasted with the holo- 
graphic Green functions computed at a finite zo, rather 
than those at the original boundary z ~ 0. 

Thus, in order to reproduce the effects of the singu- 
lar interaction (18) with £ = l,p = 2 and d — 2 in the 
holographic setting, one needs to choose the relevant pa- 
rameters as follows: r\ = 2/3 = 2(1 — a). For comparison, 
the case of £ = l,p = 1 can be covered by choosing 
r] = /3, a = 1, whereas the case of £ = 0, p = 2 requires 
X] = 1 — a, /3 = 0, which values would be unattainable 
within the class of hyperscaling violating metrics. 

Notably, all of the above metrics comply with the cri- 
teria (7) for the existence of a consistent gravity dual, 
the second one being satisfied as a strict equality, rj — 
1 - a + 13 or C = 1 + 6/d = (2(1 - a) + /3)/(l -a + 0). 
Also, the corresponding values of the dynamical expo- 
nent (£ = 3/2,1, and 2, respectively) agree with those 
inferred from contrasting the quasiparticle dispersion vq 
against the fermion self-energy^— 



E(w) 



ded d q D(u, q) 
(2ir) d+1 iuj + ie- vq 



d-l + j 



(24) 



which comparison yields £ = (£ + p)/(d — 1 + £). 

Obviously, the matching conditions (23) are only the 
necessary ones, so they may not always guarantee that 
the entire two-parameter functional dependence of the 
action S(t, x) would be reproduced with this choice of 
parameters. Our discussion should then be viewed as 
merely suggestive of a possible holographic correspon- 
dence between the aforementioned theories, and in order 
to further strengthen the case more observables would 
need to be matched. 

However, if the conditions (23) are not met then no vi- 
able gravity duals of the aforementioned physically reve- 
lant NFL systems may be found amongst the generalized 
family of metrics (5). By this argument, one concludes 
that none of such systems can be naturally mated with 
the classical AdS-RN metric considered in Refs^. 



Coulomb interacting Dirac fermions 

The recent upsurge of interest in graphene and topo- 
logical insulators - as well the earlier advent of ID 
Coulomb metals (e.g., carbon nanotutes), gapless 2D 
high-T c superconductors, and quasiparticle properties 
of 3D superfluid He 3 - brought out the problem of 
(pseudo)relativistic Dirac fermions with isolated Fermi 
points and potentially long-ranged (due to a lack of 
screening), albeit nearly instantaneous, interactions. 

In the ID case, the asymptotic dual geometry is AdS$ 
and the corresponding (conformally invariant) boundary 
theory is that of chiral ID fermions. Its proper UV com- 
pletion is naturally achieved with the use of Eq.(3) for 
p = 0, known as the 'non-rotating BTZ black hole—, 
and the (exact) finite-temperature chiral fermion propa- 
gators read 



ttT 



) 2A+ (- 



sinh nx+T 



n2A_ 



sinh(a;_T) ' 



(25) 



where x± — x ± r. In the T — > limit Eq.(26) amounts 
to G±(t, a;) = 1/x, + x_ ~ . 

According to the holographic principle^, the boundary 
theory then must be strongly coupled, as manifested by 
the UV (left /right) fermion dimensions, A± = mL/2 + 
1/2 ± 1/4, which are necessarily large, A + + A_ > 1—. 

Notably, such dimensions can not be obtained from any 
ID theory with short-ranged repulsive couplings where 
the corresponding Luttinger parameter would be limited 
to the interval 1/2 < K < 1, thereby resulting in 1/2 < 
A++A_ = i(i^ + l/if)<5/8. 

In fact, the still lower K values, < K < 1/2, can only 
be attained in the presence of long-ranged interactions, 
such as Coulomb, which endows the Luttinger parame- 
ter and the fermion dispersion with a slow momentum 
dependence: 



K(q) = 



1 



\Jl + oj Inq 



-., e q = q\/l + a\ In 



<1\ 



(26) 



where a = 2e 2 /ir. 

In the general case of the d-dimensional Dirac fermions 
with the 3D Coulomb interaction, U q ~ the coun- 

terpart of Eq.(20) reads (after the frequency integration) 



S{t, x) = J (1 - cos(e g r - qx)) 



(2nj c 9 

In the ID case, one then obtains the leading behavior 



(27) 



S(T,a?)~tr 1 / 2 ln 8 / a |i-T| 



(28) 



which gives rise to a faster-than-algebraic decay of the 
propagator G(r, x) ~ e~ s , thus implying that the system 
undergoes the ID analog of the Mott transition^. 

Making use of Eq.(15), one observes that capturing the 
behavior (28) in the holographic framework would require 
a logarithmic deformation of the asymptotic AdS% geom- 
etry 



g{z) = {L/zf In z/z Q , f(z) = h{z) - (L/z) 



(29) 
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Although this ansatz may not be unique, it shows that 
a prospective gravity dual of the Coulomb-interacting 
Dirac fermions is likely to lie outside the family of the 
AdS-RN metrics utilized in Refs<2. On the other hand, 
there exist solutions showing logarithmic behavior at in- 
termediate values of z for some generalized dilaton po- 
tentials, as in Eq.(4), tuned to their degeneracy points^. 
Thus, one would have to tap into those resources if vi- 
able candidates to the role of the gravity dual of the ID 
Coulomb metal were to be found. 
Summary 

In conclusion, we demonstrate that in order to ascer- 
tain the status of the holographic approach in the context 
of its applications to the conjectured NFLs one needs to 
venture out of the comfort zone of the customary grav- 
itational backgrounds in search of new (self-consistent) 



solutions to the coupled equations for the metric, dila- 
ton, gauge, and matter fields along the lines of Refs.— ~— . 

To that end, we discuss the specific examples of non- 
rclativistic fermions coupled to gauge-like fields and 
Dirac fermions with the Coulomb interactions, in both 
cases finding the prospective metrics to possibly belong to 
the solutions (alongside the hyperscaling- violating ones) 
of some generalized dilaton gravity models. 

Based on these observations, we suggest that within 
such a broad class of metrics one would have better 
chances of 'reverse engineering' the gravity duals of the 
already documented NFLs, thus putting the entire holo- 
graphic machinery up to a decisive test. Then, after hav- 
ing affirmed that this approach might work in the already 
known cases, one can continue expanding the list of novel 
NFL states with a greater confidence. 



1 S. A. Hartnoll, Class. Quant. G rav.26, 224002 (2009); 
J. Polchinski, larXiv: 1010.61341 J. McGreevy, Adv. 
High Energy P hys. 2010, 723105 (2010). S.Sachdev, 
larXiv:1108.1197l 

2 M. Cubrovic, J. Zaanen and K. Schalm, Science 325, 
439 (2009); larXiv:1012.568T ; S. S. Lee, Phys. Rev. D79, 
086006 (2009); H. Liu, J. McGreevy and D. Vegh, ibid, 
D83, 065029 (2011); T. Faulkner et al, ibid D83, 065029 
(2011). 

3 J. Gauntlett, J.Sonne r, D. Wa ldram, arXiv: 1106.46941 
S.Gubser, J.Ren, arXiv: 1204.63151 

4 M.Edalati et al, larXiv:1012.375T1 

5 Q. Si and F. Steglich, Science 329, 1161 (2010). 

6 S. Kachru, X. Liu and M. Mulligan, Phys. Rev. D 78, 
106005 (2008). 

7 S. A. Hartnoll and A. Tavanfar, Phys. Rev. D83, 046003 
(2011); S. A. H artnoll, D. M. Hofman, and D. Vegh, 
larXiv:1105.3197l 

8 C. Charmousis et al, JHEP 1011, 151 (2010); ibid 1201, 
089 (2012); B. Goutr aux, E. K iritsis, ibid 1112, 036, 
(2011); N. Iizuka et al. ilarXiv:1105.1162l 

9 L.Hu ijse, S.Sachdev, B .Swingle.arXiv: 1112.0573; Xi Do ng 
et al, larXiv:1201."T905l E. P erlmutter, larXiv:1205.0242l 

N. Iqbal, H. Liu, M. Mezei, larXiv:1105.462"Tl T. Faulkner 



and J. Polchinski, JHEP 1106, 012 (2011). 

1 T. Faulkner, J. Polchinski, larXivilOOl. 50491 V. Keranen, 
L. Thorlacius, I arXiv: 1204.03601 

2 S.-S. Lee, Phys. Rev. B80, 165102 (2009); M. A. Metlitski 
and S. Sachdev, ibid B82, 075127 (2010); ibid B82, 075128 
(2010). 

3 B. L. Altshuler, L. B. Ioffe, and A. J. Millis, Phys. Rev. 
B50, 14048 (1994); A. V. Chubukov, Physics 3, 70 (2010). 

4 T. Holstein, R. E. Norton, and P. Pincus, Phys.Rev.B8, 
2649 (1973); M.Reizer, ibid B39, 1602 (1989). 

5 B.I. Halperin, PA. Lee, and N. Read, Phys. Rev. B47, 
7312 (1993). 

6 J. A. Hertz, Phys. Rev. B 14, 1165 (1976); A. J. Millis, 
Phys. Rev. B 48, 7183 (1993). 

7 E. Fradkin et al, Ann. Rev. Cond. Matt. Phys.l, 153 
(2010). 

8 D. V. Khveshchenko and P. C. E. Stamp, Phys. Rev. B49, 
5227 (1994); D. V. Khveshchenko, ibid B49, 16893 (1994); 
ibid B52, 4833 (1995). 

9 D. Birmingham, I. Sachs, and S. Solodukhin, PHL 88, 
151301 (2002); D.Maity et al, Nucl.Phys. B 839, 526 
(2010); V. Balasubramanian et al. larXiv:1012.4363l 

H.Schulz, Phys.Rev.Lett.71, 1864 (1993). 



